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ABSTRACT 

We have modelled 38 barred galaxies by using near-IR and optical data from the 
^f-^ • Ohio State University Bright Spiral Galaxy Survey. We constructed the gravitational 

potentials of the galaxies from -ff-band photometry, assuming constant mass-to-light 
ratio. The halo component we chose corresponds to the so called universal rotation 

■ curve. In each case, we used the response of gaseous and stellar particle disc to rigidly 

rotating potential to determine the pattern speed. 

Q . We find that the pattern speed of the bar depends roughly on the morphological 

type. The average value of corotation resonance radius to bar radius, 71, increases 
from 1.15 ±0.25 in types SBO/a- SBab to 1.44 ±0.29 in SBb and 1.82 ±0.63 in SBbc 
. - SBc. Within the error estimates for the pattern speed and bar radius, all galaxies of 

type SBab or earlier have a fast bar (1Z ^ 1.4), whereas the bars in later type galaxies 
include both fast and slow rotators. Of 16 later type galaxies with a nominal value of 
72. > 1.4, there are five cases, where the fast rotating bar is ruled out by the adopted 
error estimates. 

We also study the correlation between the parameter 1Z and other galactic proper- 
ties. The clearest correlation is with the bar size: the slowest bars are also the shortest 
bars when compared to the galaxy size. A weaker correlation is seen with bar strength 
\Q • in a sense that slow bars tend to be weaker. These correlations leave room for a pos- 

' sibility that the determined pattern speed in many galaxies corresponds actually that 

OO , of the spiral, which rotates more slowly than the bar. No clear correlation is seen with 

■ either the galaxy luminosity or colour. 
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1 INTRODUCTION 

Studies in near-IR, where the extinction is lower than 
in visual wavelengths and majority of light comes from 
the old stellar population, have shown that about 60- 
70% of all spiral galaxies have a larg e scale stellar bar 
dEskridge et al.ll2000l: IWhvte et all |2002| ; lLaurikainen et al. " 



2004 : IMenendez-Delmestre et alj 120071 ; iMarinova fc Jogee 
20071 ). According to rece nt analysis of over 2000 spiral 



galaxies (|Sheth et alj|2008h . the bar fraction decreases from 
about 65% in the l ocal universe to abou t 20% at redshift 
0.84 (see also lAbraham et al.l Il999l; lElmegreen et all 



120041 : I Jogee et all 120041 : IMenendez-Delmestre et al.l 120071 ) 

Anyhow, bars are so common that they are either very 
robust or they represe nt a recurrent phenom enon in the 
life of a spiral galaxy (IBournaud et al.l 20051). In contrast 
with s ome earlier studies (lThompsonlll98ll : lElmegreen et al.l 
Il990l ). iHernandez-Toledo et all l|2007l ) claim the bar fre- 



quency is roughly the same in different environments - mod- 
est interactions do not seem to play a major role in bar 
formation and destruction. 

Bars can be roughly divided into two classes based 
on their light distribution: flat and exponential bars 
|Elmegreen fc Elmegreenl|l985l : lElmegreen et aLlfl996al ). In 
the first type, which is more typical to early type barred 
galaxies, the radial surface brightness profile along the bar 
major axis is natter than in the surrounding disc, whereas 
in the exponential bars the profile is quite similar to the sur- 
rounding disc. F urthermore, the flat bars can display twist- 
ing of isophotes (|Elmegreen et al1ll996bl ). Bar morphology 
can also be either classical or of ans ae-type, characterized by 
blops at the both ends of the b ar l|Laurikainen et akl 120071 : 
iMartinez-Valpuesta et al.ll2007l ). 

Another approach to characterize a bar is its strength. 
Some attempts have based on the el lipticity of the de- 
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projected bar (e.g. iMartinl 1 19951 : IWhvte et all 120021 : 
lLaurikainen fc Salo 2002]). Recently, there have been at- 
tempts to estimate the actual gravitational perturbation 



2 Rautiainen et al. 



a 




Figure 1. A schematic drawing of fast and slow bars. 



of th e bar by using near-IR p h otometry |Buta fc Blockl 
2001 ; lLaurikainen fc Said 120021 : lLaurikainen et "al ] 12004 



2005; Bu taet alj|2005j ). There is also another bar strength 
estimate, namely the ^-Fou rier amplitudes of density 
(Laurikaincn et al ] 12004 l2005h . which is an approximation 
of the relative mass of the bar. All these bar strength esti- 
mates are discussed with respect to the Hubble sequence by 
lLaurikainen et al.1 l)2007t ). 

Perhaps the most important parameter defining a bar 
is its pattern speed, Qbar, or how fast the bar rotates. In 
principle, this determines how far the orbits of stars and gas 
clouds are affected by the bar. The pattern speed has a phys- 
ical upper limit - a bar cannot reach beyond its corotation 
resonance (CR) radius Rcr, i.e. the region in the disc where 
the angular speed of circular rotation equals the bar pattern 
speed. This limitation is based on the studies of stellar or- 
bits in barred potentials - the orientation of the major axes 
of closed orbits becomes perpendicular with the bar beyond 
Rcr, thus th e orbits in the out er disc are not able to sup- 
port the bar (|Contopoulos|[l980l ) . On the other hand, there 
is no evident lower limit based on stellar orbits for the bar 
pattern speed. 

In the literature an often used nomenclature is based 
on the value of a dimensionless parameter 1Z = Rcn/Rbar, 
where Rbar is the semi major axis of the bar. The cases where 
1Z ^ 1.4 are usually called "f ast bars" , whereas those with a 
larger ratio are "slow bars" (|Debattista fc Sellwoodl fcoOO). 
This is shown as a schematic drawing in Fig. [T] Consider- 
able effort has been devoted to determine 1Z for individual 
galaxies and to study if it depends on other properties of 
the bar itself or the other galactic components. 

One way to study the evolution of galaxies is to make 
self-consistent N-body simulations, where the original parti- 
cle distribution forms different kinds of structures depending 
on its initial stage. If a bar forms by a global bar i nstability 
in th ese simulations, it tends to be fast rotating (|Sellwoodl 
Il98ll ). On the other hand, it has been suggested that if a bar 
forms in a galaxy due to an interaction with other galax- 
ies, it can be a slow rotator, perhaps extending onl y near 
its inner Lindblad resonance (tMiwa fc Noguchlll998l). This 
may also be the case with gradual bar growth (|Lvnden-Belil 



or when the initial b ulge-to-disc mass ratio is low 
1 Combes fc Elmegreenlll993l ). The pattern speed of the bar 
does not depend only on its initial value - the bar interacts 
with the other components of the galaxy. Secular evolution 
due to interaction between the bar and the outer disc decel- 
erates the bar rotation, but this seems to be compensated 
by the growth o f the bar length by capturing more particles 
jSellwoodlll98lh . The situation is different with the interac- 
tion between the bar and the dark halo: if the halo density 
is high in the inner parts of the galaxy, the resonant inter- 
action can decelerate an initia lly fast rotating bar so much 
that it becomes a sl ow rotator (|Debattista fc Sellwoodiri998l ; 
lAthanassoulal l2003). Thus, the observed bar pattern speed 
could be considered as a rough estimator for the halo con- 
tribution. 

Much weight has been given to a model-independent 
kinematic method, the so-called Tremaine-Wein berg 
method, hereafter TW (|Tremaine fc Weinberg! I1984T I. It 
is based on photometric and spectroscopic observations 
along slits parallel to the galaxy's major axis and assumes 
a rigidly rotating component, which satisfies the continuity 
equation. TW-method has been used to determine pattern 
speeds of about 20 large scale stellar bars. The results range 
from probably unphysical 1Z = 0.6 to slow bar regime with 
1Z = 1.8. This has been considered to be in accordance 
with fast bars, at least when the quite large error estimates 
(e.g. due to sensitivity to errors in the position angle of the 
disc. iDeba ttista 2003 ) are t aken into account (lKentlll987l; 



i Merrifield fc Kuiikenl 1 19951 ; 
19991 . 120031; 
20031: 



Gerssen, Kuijkcn & Mcrrificld 



Aguerri. Debattista fc Corsini 
"20041 : 



Debattista & Williams 
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20071 : iTreuthardt. Buta. Salo fc Laurikainenl I2007T I. This 
apparent lack of slowly rotating bars has been interpreted 
to indicate that the density of a dark halo is small in 
the inner parts of barred galaxies (|Debattista fc Sellwoodl 
l2000l : lAguerri et al] 120031 ), which is in disagreement with 
some cosmological simulations of structure formation, typ- 
ically producing cen trally peaked dark matter component 
(jNavarro et al.lll996t ). However, due to limitations set by the 
initial assumptions, measurements using stellar absorption 
line spectrometry have been made almost exclusively for 
SBO-galaxies, which are almost free of gas and dust. 

Although som e attempts l|Zimmer et al.l 120041 : 
iRand fc Wallinl I2004T ) have been made to use TW-method 
with CO observations, the pattern speeds of late type 
barred galaxies have been usually determined by various 
indirect methods. Many of these rely on morphological 
features in barred galaxies such as spiral arms, rings or 
leading offset dust lanes. In a classic scenario, a two-armed 
spiral starts from the ends of the bar, but exceptions to 
this are not rare: there can be an "empty" region between 
the bar and the spiral (e.g. NGC 210) or an offset between 
the bar position angle and the starting points of the spiral 
arms (e.g. NGC 799). There can also be more than two 
arms either in the whole scale of spiral structure (e.g. ESO 
566-24) or as a multiply armed outer disc (e.g. NGC 4303). 
The clearest correlation with the bar properties is observed 
in t he sizes and shapes o f rings (for a detailed review, 
SCO iButa fc Combes! Il996l ): outer rings whose semi major 
axis is about twice the bar radius and inner rings whose 
major axis is usually the same as the major axis of the bar. 
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There is also a third ring type, the nuclear rings, whose 
radius is roughly 1/10 of the bar radius, but with a large 
scatter. Leading offset dust lanes are straight or curved 
dust features inside the bar and are located in its leading 
side. 

Models of individual galaxies have been constructed ei- 
ther by fitting analytical potential components to observa- 
tions or by determining the potential from photometry, mak- 
ing some assumption about mass-to-luminosity ratio, orien- 
tation and internal geometry of galaxies. In principle, mod- 
elling is more economical than using TW-method - the pat- 
tern speed can be estimated by comparing the simulated and 
observed morphology. Naturally, kinematical observations 
can help to determine the acceptable parameter range. Mod- 
elling has produce d both high and low pattern speed esti- 
mates for the bars (iHunter et al.lll988l; ISempere et al. 19951; 
Lindblad fc Kristenlll996l : ISalo et aLlll999l : iRautiainen et all 
20041 ). The given error estimates of models are typically 
smaller than with TW-method. 

Pattern speed estimates have also been done by iden- 
tifying various morphological or photometric features with 
resonances. Inner and outer rings are usually considered 
to form by gravitational torque of the rotating bar, which 
causes gas to flow in the radial direction. The net torque 

vanishes in major resonances, where gas then accumulates 
■ " — it — ■ — | 

(jButa fc CombesHl996r ). Comparison with results of gas dy- 
namical simulations and analysis of orbits in barred galax- 
ies have led to the following identifications: an outer ring 
should be located near the Outer Lindblad resonance (OLR) 
and the inner ring nea r the inner 4/ 1 -resonance, which 
is located inside Rc r dSchward Il98ll. Il984l : iBvrd et al ' 



19941 ; ISalo et al.lfl999l : IRautiainen fc Saloll2000l ; iPatsis et al 



2003). When the rotation curve is known, these resonance 
identifications can be used to determine the pattern speed 
l|Buta fc PurcellHl998fh Another approach is based on two- 
colour photometry (jPuerari fc Dottoril 1 19971 ) : the location 
of shock induced star forming regions should change side re- 
lated to spiral arms when crossing Rcr- This crossing should 
be detect able from two - colour photometry. Especially, the 
results of lAeuerri et all (|l998h gave a hint (the total sam- 
ple consisted of 10 galaxies) that late type spirals could be 
slower rotators than early types. 

In addition there are at least two other methods worth 
mentioning. A kinematical method to determine the corota- 
tion radius, which is based on the residual patterns in the ve- 
lo city field after re moval of circular velocities, was suggested 
bvlCanzian (19 9jl). However , it has been rarely applied. Re- 
cently, IZhang fc Butal I^OOTf ) suggested that the calculated 
phase shifts between potential and density could be used to 
determine Rcr- The validity of this approach is still ques- 
tionable: for several galaxies it found corotation radii well 
inside the bar, which is in disagreement with the analysis of 
orbits in barred potentials. 

Using the previously published pattern speed estimates 
in studying a possible dependency on the morphological type 
is problematic. Very different methods have been used, the 
data is of uneven quality or the definition of the bar differs. 
There are several papers where quite successful simulation 
models are presented, but only one pattern speed has been 
tried (e.g. by assuming TL = 1) - this cannot be taken as a 
determination of the pattern speed. 

Large galaxy surveys and the increase of computing 



power makes it possible to improve the situation by mass 
production of galaxy models: here we present a simulation 
series where we estimate the pattern speeds of 38 moderately 
inclined barred galaxies, using data fro m the Ohio State Uni- 
versity Bright Spiral Galaxy Survey jEskridge et all |2002| . 
hereafter OSUBSGS). To our knowledge, this is the largest 
sample of barred spiral galaxies whose pattern speeds are de- 
termined with a consistent method. Furthermore, the mor- 
phological types of these galaxies range from SBO/a to SBc, 
based on Third Reference Catalogue of Bright Galaxies 
|de Vaucouleurs et alj [l99ll . hereafter RC3), so that each 
morphological type is represented by several galaxies. The 
initial modelling r esults for this sample were published in 
IRautiainen et all (|2005l . hereafter RSL2005), a nd part of 
the an alysis in this paper was also presented in ISalo et al.l 
(|2007h . 

Whereas we find all the bars of early type galaxies 
(SBab or earlier) of the sample are fast rotators, our mod- 
els also suggest that galaxies of later morphological types 
include both fast and slow rotating bars. The slow rotation 
seems to be related to the small size of the bar. An alter- 
native interpretation to this is that the spiral arms in these 
galaxies have a lower pattern spe ed than the bar, a situa- 
tion often seen in N-bod y models ijSellwood fc Sparkdll988l : 
IRautiainen fc Salolll999l ). 



2 MODELLING 

Our modelling method is essentially the same as we used 
in modelling IC 4214 and ESO 566-24 (|Salo et all Il999l : 
IRautiainen et alJliooi l. the main difference being that we 
do not have kinematical data for this larger sample. We as- 
sume that the //-band light distribution corresponds to the 
projected distribution of luminous matter (constant M/L) 
and that there is only one pattern speed in each galaxy. 

We simulated the behavior of two-dimensional discs of 
collisionless stellar test particles and inelastically colliding 
gas particles in the determined potentials. In each collision, 
the relative velocity component of two gas particles in the 
direction joining the particle centres was reversed and mul- 
tiplied by the coefficient of restitution, here assumed to be 
zero. The main modelling parameter was the pattern speed 
of the non-axisymmetric part of the potenti al. Fo r more de- 
tails on the simulation code, see ISalol (|l99lh and lSalo et ail 
( 1999). The best fitting pattern speed and a rough estimate 
of the error (the average error estimate for Rcr is about 
±15%) was determined by visual comparison between the 
simulation morphology and the morphology in the B- and 
Jf-band images. 



2.1 Mass model 

We used the //-band images from the OSUBSGS, for 
which a bulge-disc-ba r decomposition was made by 
lLaurikainen et all (2004 S). This multiple-component process 
gives more realistic bulge parameters than fitting just disc 
and bulge. For each image the bulge component was first re- 
moved, after which the disc was deprojected to face-on ori- 
entation. The light distribution and surface density of the 
disc was approximat ed by a Fourier decomposit ion using the 
same approach as in lLaurikainen fc Said l|2002l ). 
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Figure 2. The calculated rotation curve (top) and corresponding 
frequency diagram (bottom) for NGC 4303. The continuous line 
in the top frame shows the rotation curve, the dotted line the 
bulge contribution, the dashed line the disc and dot-dashed line 
the halo contribution. In the bottom frame the curves drawn with 
continuous line corresponds to Q and Q ± re/2, the dashed lines 
corresponds to Q ± re/4 and the the dotted line indicates the best 
fitting pattern speed. Also the bar radius Rbar ~ 53" is indicated. 
This close coincides with inner 4/1-resonance. 



The disc gravity was calculated using even azimuthal 
components from m=0 to 8. The strongest of the odd com- 
ponents is usually m=l. However, it can represent lopsided- 
ness caused by a recent interaction and even if it is related to 
an internal mode, it can have a different pattern speed than 
the xa=2 mode. The omission of odd modes can cause some 
of the models look more regular than the observed galaxy. 

In the vertical direction, an exponential distribution 
was assumed, with a constant scale height throughout the 
disc: the scale height was chosen according to the morpho- 
logical type of the galaxy, using the typi cal ratio of ver - 
tical to radial scale lengths ha given bv Ide Griisl (fl99^l . 
The gravitational potential of the bulge was added to the 
disc potential, assuming that the bulge mass is spheri- 
cally distributed. In our standard mass model we also in- 
cluded a dark halo component ba sed on the universal ro- 
ta tion curve of iPersic et alj |l996), in a similar manner as 
m iButa et all l|2004 l. To calculate the halo profile, we used 
L/Li, calculated from the RC3 and NED (NA SA/IPAC Ex - 
tragalactic Database) data and distance from iTullvl (Il988l ) 
(except for NGC 6782 for which the distance was calcu- 
lated from the radial velocity assuming Hubble constant 
H = 75 km s" 1 Mpc" 1 )- For more details of the poten- 
tial calculation see iLaurikainen fc Said l|2002ft . An example 
of a calculated rotation curve and corresponding frequency 
diagram showing £1 (circular frequency) and Q±k/2 (where 
k is the epicycle frequency) , is shown in Fig. [2] 

Because we do not have kinematical data, we cannot 
fine-tune our mass models. However, we can estimate the 
robustness of our pattern speed fits by trying different mass 
models. In addition to our standard model with universal 



rotation curve, we also made simulations without the dark 
halo component. In general, the calculated rotation curves 
were then falling in the outer parts of the discs, but the 
determined pattern speeds were quite similar because the 
fitting was done mostly using features in the inner parts. 
For eight galaxies (six with slow bars, two with fast bars), 
we also studied the effect of increasing the halo contribution. 
This was tested by multiplying the even Fourier amplitudes 
m = 2 — 8 of the force by a factor 0.75 or 0.5. The effect to 
the fitted pattern speed was quite small, in all cases within 
the error estimates of the standard model series. In a couple 
of cases, the fit was marginally better when the amplitudes 
were multiplied by 0.75, but in all cases the fit was worse 
with multiplication factor 0.5. This is in good agreement 
with our tests wit h the mass models for ESQ 566-24 (see 
Figures 4 and 12 in iRautiainen et alll2004h . where we found 
that three mass models (no halo, a halo closely correspond- 
ing to universal rotation curve and dominating halo) all gave 
essentially the same pattern speed, although the quality of 
the fit was clearly different. Also, modifying the mass model 
so that the major axis rotation curve determined from the 
simulated velocity field matched the observed one as closely 
as possible, changed the corotation resonance distance only 
slightly. 



2.2 Bar radius 

There is no unambiguous way t o determine the bar radius 
(for a detailed analysis see e.g. lAthanassoula fc Misiriotisl 
[2002J). Some approaches rely on Fourier analysis of the sur- 
face brightness: the bar is identified with the region where 
the m = 2 (or m = 4) phase angle is almost constant, or 
where the amplitude drops to a certain fraction of its maxi- 
mum value. Another method is to fit ellipses to isophotes - 
bar radius is then determined as the radius where the ellip- 
ticity reaches a maximum or where there is a steep drop in 
it. Besides these, there are also several other methods. 

In RSL2005 we estimated lengths of the bars by us- 
ing visual inspection of deprojected H-band images and 
isophotes and m = 2 Fourier phase angles. In this st u dy, we 
have basically followed the same approach as Erwinl (|2004| . 
12005ft (see also lMichel-Dansac fc Wozniakll2006ft , with slight 
modifications. We first remove the surface brightness of the 
fitted bulges from the images and then deproject them. Next 
ellipses are fitted on the isophotes of the deprojected images 
by using routine ELLIPSE in IRAF. For the lower limit of 
the bar radius we normally use the semimajor axis of max- 
imum ellipticity a e (the outermost maximum that is still 
within the bar). As the upper limit we use Lj, ar which is the 
minimum of two ellipse-fit measures - the first minimum el- 
lipticity outside a e , or the point where the position angles 
of the fitted ellipses differ more than 10° from the position 
angle of the bar (the smaller one of these is chosen). If a 
galaxy have an inner ring that interferes the determination 
of a E and/or L bar , we use the inner edge of the ring as the 
lower limit and its outer edge as the upper limit for the bar 
radius, both measured along the bar major axis. The bar 
radius is then determined as Rbar = (a e + Lf )ar )/2 and its 
error estimate is taken as ARtar = Rbar — a e = Lbar — Rbar- 
The adopted values are given in Table 1, which also lists the 
error estimates. The error estimates in the bar radius are in 
range ±4% - ±24 %, the average being about ±14%. 
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When comparing our results with other bar length mea- 
surements, we limit to those done with near-IR photom- 
etry. The bar lengths used in this stu dy are on average 
about 10% longer than the values given bv lLaurikainen et al.l 
|2004 ). who used the phase angle method. This typically 
corresponds to our lower limits for the bar lengths. In five 
ca ses our bar radius estim ate is clearly smaller compared 
to lLaurikainen et al.l (|2005h : the difference is related to the 
identification of the bar: for example if a galaxy has an outer 
ring of type i?i , then the innermost part of the ring is close 
to the major axis of the bar, and can make the m = 2 phase 
angle almost constant far beyond the actual bar region. A 
good example of this is NGC 3504. 

Eleven of our galaxies are common with lErwinl l|2004 
l2005h , who used essentially the same method to determine 
the bar radius. Their bar radius estimates (taking into ac- 
count possible differences in the adopted orientation param- 
eters) are in good accordance with our results: typically our 
bar radius estimate is between their values for a e and Lbar, 
and if that is not the case, the difference is only marginal. 
The only exception of this is NGC 4303, for which we get 
60 % higher value. 

IMarinova fc Jogeej (|2007f t used the OSUBSGS images 
in their study of bar properties. They give bar radii based 
on isophote ellipse fitting for 34 galaxies that are also in our 
sample. Typically our bar radius estimates are few percent 
larger than their values (median is about 5%), which is not 
surprising taking into account that this is expected to corre- 
spond to ou r lower limit. The differen ce is largest for NGC 
1317 where IMarinova fc Jogeel (|2007l ) gives measurements 
for a feature we consider to be a secondary bar or a nuclear 
ring. It is the only case in their sample, where the given 
bar radius is not within our error estimates. In addition to 
previous ex a mples , three of the galaxies are common with 
lLaine et al.l l|2002T ) (who used the ma ximum of fitted ellip- 
ticity) and one with lReese et al.l (|2007t ) (distance where 10% 
of the maximum of fitted bar intensity profile is reached). 
Our bar radius estimates agree well with their results when 
error estimates are taken into account. 



2.3 Defining the best fit 

The morphological features we used in the comparison are 
typically the spiral arms, the inner rings which surround the 
bar, and the outer rings. With spiral arms we compared their 
extent, location and pitch angles, and occasionally other fea- 
tures such as straight arm segments. With rings we com- 
pared the shapes and sizes to the observed ones. This com- 
parison was done with the help of a crosshair figure (cross 
+ circles indicating the disc scale length and its multiples), 
that was overplotted on both the galaxy images and the 
modelled gas and stellar particle distributions. 

An example of our modelling approach is shown in 
Fig. |31 where we compare models with six different pattern 
speeds to H- and B-band images of NGC 4303. In a few 
cases we could not determine a single best-fitting pattern 
speed because the outer parts favored lower pattern speed 
than the inner parts. A possible reason for this is that a 
galaxy might have many pattern speeds so that the spiral ro- 
tates more slowly than the bar. T his is often seen in N-bod 
simulations of barred galaxies (ISellwood fc Sparkd Il98 
iMasset fc Taggerlll997l ; iRautiainen fc Salolll999r ). If two fits 



are considered to be of equal quality, the value adopted 
in this study corresponds to the smaller TZ (faster bar), 
best characterizing the bar region. The ambiguity of pattern 
speed determination was also faced bv lLindblad fc Kristenl 
(1996): for galaxy NGC 1300 they found two values giving 
TZ ~ 1.3 or TZ « 2.4, both producing equally good fit. 

Several of the galaxies contain leading offset dust lanes 
in their bars. Although such features are sometimes seen 
also in our models (see Fig. [3)l , we did not include them to 
pattern speed determination: the dust lanes correspond to 
highest density in the interstellar medium an d it is probable 
that our gas modelling method (|Sald ll991) is not suitable 
for such situations. 

Most of our models contain features in their gas com- 
ponents that look like nuclear rings observed in galaxies. 
Although several studies have connected nuclear rings to 
inner Lindblad resonance (ILR) or resonances, it has been 
claimed, based by gas dynamical models, that they are ac- 
tually not related to resonances but to orb its in the cen- 
tral part of the bar (|Regan fc Teubenll2003r i. Furthermore, 
in the models of Regan & Teuben the nuclear rings shrank 
even when the gravitational potential was constant. Even 
so, in our previous models of IC 4214 and ESO 566-24, the 
model which gave the best fit to overall morphology and 
kinematics, produced also a nuclear ring whose major axis 
diameter was the same with that of the observed ring. How- 
ever, we have not included the nuclear rings to comparison 
with the sample galaxies: the image quality is usually not 
good enough to show a possible nuclear ring and the grav- 
itational potential, calculated from the //-band image, is 
probably not accurate enough in the innermost parts due 
to low resolution (one pixel corresponds to 0.05 - 0.2 kpc) 
Also, some of the galaxies included active galactic nuclei 
or central star bursts, which affect the central parts of the 
brightness profiles even in the //-band, and thus also the 
deduced potentials (see Figs. 2 and 3 in lSalo et al.|[T9 99). 

In addition to gaseous nuclear rings, our models of- 
ten show distinct stellar features, whose major axis is 
perpendicular to the bar, in the same region. These fea- 
tures, which are us ually n o t seen in real galaxies, are re- 
lated to X2-orbits. Patsis (2005) demonstrated that X2- 
features can be removed from the collisionless test par- 
ticle models by fine-tuning the initial conditions. On the 
other hand, this region often has a secondary bar of about 
similar size as a typical nuclear ring b ut a different pat- 



tern speed than the main bar (e.g. Frie dli & M artinet] 



19931: IRautiainen fc Salolll999l; iMacieiewski fc Sparkd |2000|: 



Heller et al.l l200ll : IRautiainen et all l2002i : iMacieiewskil 
20061 ). 



3 NOTES ON INDIVIDUAL GALAXIES 

Here we give comments on some of the modelled galaxies. 

NGC 289: this SBbc (T = 4) galaxy has a small bar, 
the inner spiral arms start from the ends of the bar. There 
is also an extended outer spiral structure, which bifurcates 
from the inner spiral arms. The best-fitting model, which 
sets TZ « 2.6 reproduces the inner spiral, the bifurcations 
and the inner part of the outer spiral (see Fig. |4}. The out- 
ermost part of the outer spiral structure is not reproduced, 
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Figure 3. Models with different pattern speed for NGC 4303. First two frames shows deprojected images of the galaxy in H- and 
B-bands. The number in the simulation frames indicates TZ = Rcn/Rbar- The radii of the drawn circles corresponds to 1, 2 and 3 scale 
lengths of the disc in the //-band. 



but it can be related to an interaction with a small compan- 
ion galaxy, LSBG F41 1-024. model sets Rcr much outside 

NGC 613: this SBbc (T = 4) galaxy has a bar, inner 
ring and several spiral arms emerging from the ring. The 
Fourier phase angle method gives too high bar radius, well 
in the area of spiral arms. The fit, which sets TZ ~ 1.6, is not 
very good (see Fig. [4| : only the inner ring and small parts of 
the arms are reproduced. The inner ring would allow also a 
bit higher pattern speed, which would set the galaxy clearly 
to the domain of fast bars. 

NGC 1187: this SBc (T = 5) galaxy has a bar with 
three spiral arms emerging from it. The best-fitting model 
(see Fig. [4}, setting TZ ~ 2, has four arms due to using 
even Fourier components, three of which corresponds well 
the observed arms. The error estimate for TZ is quite large, 
and a fast-bar solution cannot be ruled out. 

NGC 1241: in this SBb (T = 3) galaxy the bar is 
surrounded by an inner pseudoring, which exhibit rectangu- 
lar features, the spiral structure is not regular, possibly due 
to interaction with NGC 1242. The best-fitting model (see 
Fig. [4]), setting TZ » 1.4, reproduces well the rectangular 
structure of the inner pseudoring and also part of the spiral 
structure. 

NGC 1302: the only distinct feature in this SBO/a 
(T = 0) galaxy is the very tightly wound spiral structure. 
The best-fitting model (see Fig. [4j| , setting TZ ~ 1.6, is cho- 
sen by the extent of the gas spiral. 

NGC 1832: this SBbc (T = 4) galaxy has an inner 
ring whose major axis seems to be larger than the bar. The 
spiral structure asymmetric. The ring and one spiral arm 
are reproduced by the best-fitting model (see Fig. [4}, setting 
TZ » 1.7. 



NGC 3275: this is a SBab (T = 2) galaxy with an inner 
ring and asymmetric spiral structure. The ring and part of 
the spiral structure is reproduced with TZ ta 1.5 (see Fig. Q. 
The Fourier phase angle method clearly overestimates the 
bar radius: it is outside the lopsided inner ring. 

NGC 3504: this is a SBab (T = 2) galaxy with an 
outer ring. The ring and the extent of the bar is best re- 
produced with a simulation of stellar test particles giving 
TZ ~ 1.2 (see Fig. 3}. In the gas component the same pattern 
speed is favored, but the fit is worse because the outer ring 
does not form. The outer pseudoring in this galaxy causes 
the Fourier phase angle criterion to work incorrectly, giving 
too high value for bar length. 

NGC 3507: this SBb (T = 3) galaxy has really long 
spiral arms that can be followed about 360° . The best-fitting 
model, setting TZ ~ 1.2, reproduces the spiral structure in 
its full extent. Overall, this is one of the best models (see 

Fig. ed . 

NGC 3583: this SBb (T = 3)galaxy has an inner pseu- 
doring with rectangular features and two-armed spiral. The 
best-fitting model, with TZ m 1.2, reproduces the shape of 
the pseudoring and also the spiral structure is well fitted 
(see Fig. 0}- 

NGC 3726: this SBc (T = 5) galaxy has a weak and 
small bar. The inner part of the spiral structure is somewhat 
polygonal. The best-fitting model reproduces this appear- 
ance and gives TZ~2 (see Fig. [4|. 

NGC 4303: this SBbc (T = 4) galaxy has a small 
bar inside far extending spiral structure whose inner part 
shows straight arm segments, sometimes called "Vorontsov- 
Velyaminov rows" l|Chernin et al.|[2000l ). This structure is 
reproduced with a narrow pattern speed range (see Fig. [3j 
by a model giving TZ ~ 1.7. 
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Figure 4. Dcprojected image of the sample galaxies (B-band) vs. models (gas or stellar component). The Hubble stage T is also given. 




Figure 4 - continued 
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Figure 4 - continued 
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NGC 4548: this SBb (T = 3) galaxy has spiral arms 
that seem to start offset to the bar in the leading side - at 
least the part connecting spiral to the ends of the bar is very 
weak. The best-fitting model gives TZ « 1.3 (see Fig. [3J . One 
arm of this galaxy shows straight segments, unlike the case 
of NGC 4303, these are not reproduced by the model. 

NGC 5701: the outer pseudoring of this SBO/a (T = 0) 
galaxy is almost detached from the bar. It is reproduced with 
a model giving 72. w 1.4 (see Fig. 3}. 

NGC 5850: the inner ring of this SBb (T = 3) galaxy is 
reproduced with a model giving TZ « 1.4 or Rcr ~ 105±12" 
(see Fig. [4} . The outer structure of this galaxy is not regular, 



possi bly due to interaction with NGC 5846 (|Higdon et al.l 
1998). A hydr odynamical model for this galaxy has been 
published by lAguerri et all (|200lh . giving R CR « 90". We 
found that an alternative model, giving Rcr ~ 82 ± 6" or 
TZ ~ 1.1 is also possible when matching the size of the inner 
ring, but then the shape is not correct. 

NGC 6384: this SBbc (T = 4) galaxy has a small bar 
and long extending spiral structure that winds about 360°. 
The best-fitting model, giving TZ « 2.3 reproduces the spiral 
in its full length (see Fig. [4|. 

NGC 6782: the inner ring and the large part of the 
spiral of this SBa (T = 1) galaxy are reproduced by a model 
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giving TZ « 1.3 (see Fig. [4}. The main difference between the 
observed and modeled morphology is that the real galaxy 
has an R[ outer pseudoring whereas the spiral in the model 
does not form the pseudoring. The Fourier phase method 
clearly overestimates the bar radius by setting it well out- 
side the inner ring. This is possibly caused by the outer 
pseudoring that affects the phase angle of the Fourier de- 
composition. 

NGC 7723: the inner pseudoring and multiple-armed 
spiral of this SBb (T = 3) galaxy are best reproduced by 
a model giving TZ « 1.3. There is an alternative model, 
with TZ w 1.9, which gives about as good fit with the spiral 
structure. For comparison see Fig. [5] 



4 RESULTS AND DISCUSSION 

The results - the bar radii and the fitted corotation reso- 
nance distances are listed in Table 1, along with error esti- 
mates for both parameters. The differences to preliminary 
values for Rcr that were given in RSL2005 are minor, within 
the error estimates. With some galaxies the values of TZ has 
changed considerably due to different way to determine the 
bar radius. However, in the sample level this does not in- 
troduce any significant changes in the observed dependency 
e.g. on the morphological type. 

In addition to TZ vs. Hubble stage comparison that was 
already done for the preliminary values in RSL2005, here we 
search possible correlations between the determined pattern 
speed and the size and the strength of the bar. We also study 
if the fitted TZ depends on the bulge contribution, galaxy 
luminosity or colour. 



4.1 Comparison with previous pattern speed 
estimates 

Before going to analysing the results of the sample we first 
compare our pattern speed determinations with previous es- 
timates for these galaxies. From literature we found pub- 
lished pattern speed estimates for eight galaxies in our sam- 
ple. Unfortunately, none of these was based on direct TW- 
method, instead either morphological arguments or simula- 
tion modelling was used to determine 1Z (for NGC 1073 and 
NGC 4123 we found both morphological and model based 
estimates of pattern speed). 

In three cases, NGC 1073, NGC 41 23, and NGC 5921 
Rcr w as estimated by the method of IPuerari fc Dottorl 
1 19971 ) l|Aguerri et all Il998h : the B- and /-band m = 2 
phase angle crossing was taken to be signature of CR, giving 
1Z = 1.17, 1.41 and 1.28, respectively. These compare quite 
well with our estimates (72 = 0.95 ± 0.24, 1.17 ± 0.28 and 
1.25 ±0.23). 

For NGC 613 dEfmegreen et al.lll992h and NGC 3504 
l|Kennev et al.l Il993 ) the morphological signatures of the 
outer Lindblad resonance (extent of spiral structure, outer 
ring) and rotation curves were used to determine Rcr (in 
both cases correspo nding to 72 = 1.0). For NGC 7723 
lAguerri et al.l l|2000l ) determined 72. = 1.0 by assuming that 
the corotation region lack recent star formation. Of these, 
the cases of NGC 3504 and NGC 7723 are in good agreement 
with our value (1Z = 1.19 ± 0.27, 72 = 1.31 ± 0.29), but for 
NGC 613 the difference is large (our value TZ = 1.61 ± 0.24). 
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Figure 6. The best fitting TZ vs. Hubble stage T (decimal num- 
ber values from Hyperleda are used in this and other figures to 
increase the clarity of the plot by spreading the distribution data 
points in horizontal direction). Dashed line shows the upper limit 
for fast bars. 



However, our model of NGC 613 was not particularly good: 
the fit is based just on the inner pseudoring and the inner- 
most part of the spiral arms. Perhaps the complicated spiral 
structure in beyond our simple model. 

Three galaxies of the sample hav e been previously mod- 
elled. For NGC 1073 lEngland et ail (|l990h fo und TZ = 1.1 
four v alue TZ = 0.95 ± 0.24), for NGC 4123 IWeiner et all 
|200ll ) fou nd TZ = 1.35 (our val ue TZ = 1.18 ± 0.26) and for 
NGC 5850 lAguerri et all (|200lT l found TZ = 1.35 (our value 
TZ = 1.39 ± 0.29). In addition to these, the pattern speed of 
NGC 4314 was determined by matching the shapes of parti - 
cle orbits with the observed outer ring l|Quillen et al.lll994l ) 
giving TZ = 1.17 (our value TZ = 0.99 ± 0.19). 

Altogether, in half of the eight cases a smaller value of 
TZ was found than in our study, and in another half a larger 
one was reported in the literature. For only one galaxy the 
difference was clearly larger than our error estimate for 72. 
If we take into account differences in adopted bar lengths 
and orientation parameters, our models can be considered 
to be consistent with previous pattern speed estimates for 
these galaxies. 



4.2 Dependence on Hubble type 

Fig. [6] shows the determined values of TZ as a function 
of H ubble stage T, taken from HyperLeda l|Paturel et alj 
2003). There seems to be a trend between the pattern speed 
and the Hubble stage T: the average value of TZ and its stan- 
dard deviation are 1.15±0.25 for SBO/ato SBab (T = 0-2), 
1.44 ±0.29 for SBb (T = 3) and 1.82 ±0.63 for SBbc to SBc 
(T = 4 — 5). Thus, early type (SBab or earlier) galaxies 
are typically fast rotators (in three cases the nomimal value 
TZ > 1.4, but the fast bar is within the error estimates), 
whereas later type galaxies include both fast (10 cases) and 
slow rotators (16 cases). Especially, with five galaxies the 
fast rotating bar is ruled out based on the adopted error 
estimates for the corotation resonance radius and bar size. 
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Figure 5. Models with different pattern speed for NGC 7723. First two frames shows deprojected images of the galaxy in H- and 
B-bands. The number in the simulation frames indicates 7Z. The radii of the drawn circles corresponds to 1. 2 and 3 scale lengths of the 
disc in the _ff-band. 



Our results can be compared with other pattern speed 
estimates also in the sample level - did we find a similar 
trend as the other researchers. For this purpose we have 
composed a comparison sample of pattern speed determi- 
nations from literature by the following criteria: 1) using 
the TW-method, 2) using simulation modelling, 3) using 
the Puerari-Dottori -method, and 4) galaxies whose pattern 
speed estimate is based on morphological arguments. Each 
galaxy enters the sample only once: if its pattern speed is 
estimated by TW-method, then its pattern speed estimates 
by other methods are omitted and so on. If essentially the 
same method is used several times for a galaxy, then the 
newest measurement is selected. The comparison sample is 
listed in Table 2 and the comparison with our modelling is 
illustrated in Fig. [7] 

Strictly speaking, our results are in agreement with 
those obtained with TW-method. Most of the measurements 
with TW-method have been limited to early type galax- 
ies, which indeed have fast bars. This has been taken to 
indicate that the dark halo contribution in these galax- 
ies i s so low that it cannot d e celerate the bar r ota- 
tion l|Debattista fc Sellwoodl l200d : lAguerri et aT]|2003l ). If 
so, then the slow bars in our sample could be systems 
where angular momentum transfer between the bar and 
the dark halo with a substantial central density has slowed 
down bar rotation. Alternatively, the difference in pattern 
speeds can be primordial: the early type galaxies, which 
have massive bulges, may fa vor the formation of fast bars 
|Combes fc Elmegreenll 19931 ). 

In future it would be important to increase the overlap 
of pattern speed estimates based on different methods both 
by extending the current modelling to earlier type galaxies 
(e.g. by using the data of the ongoing survey of SO-galaxies, 



lLaurikainen et al. 2005; Buta et al1l2006l ) and by extending 
the direct measurements to later types. This would give a 
better understanding of accuracy and possible bias of pat- 
tern speed determinations done with different methods. 

The slow bars in our sample cannot be refuted by the 
ambiguity in the determination of bar length: our values 
tend to be slightly larger than previous bar length estimates 
for the same galaxies, thus favoring a faster bar. Further- 
more, the comparison with other estimates of corotation ra- 
dius for the sample galaxies does not show any systematic 
difference. If our error estimates for Rcr and Rtar (as listed 
in Table 1) hold, then of 19 galaxies with nominal value 
7Z > 1.4, we can coin five (six if we adopt the alternative 
pattern speed solution for NGC 7723) as "definitely slow ro- 
tating", i.e. 1Z — A7Z > 1.4 (see Fig. [BJ, and several others 
for which a fast bar solution in unlikely. Of these, NGC 289 
and NGC 578 both have a very small bar when compared to 
the extent of the spiral structure. The innermost structure 
of these galaxies would allow also higher pattern speed than 
the best fitting model, but even then 1Z would be much larger 
than 1.4. In NGC 1832 the best-fitting model, which gives 
a slow bar, reproduces the inner ring. The ring disappears 
when the pattern speed is increased to fast bar domain. The 
bar does not fill the inner ring, thus also the morphology 
supports slow bar solution. The innermost structure of NGC 
3726 allows also a faster bar pattern speed, but even the in- 
ner ring morphology is better reproduced with the value 
listed in Table 1, which gives a good overall fit. The straight 
arm segments of the inner spiral structure in NGC 4303 are 
well reproduced with the listed value of 7Z, but disappear 
when reaching the fast bar domain. However, the bar of this 
galaxy has leading offset dust lanes, which are usually taken 
as a signature of fast-rotating bar. Our gas dynamical model 
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Table 1. The measured bar length R^ ar and its error estimate AR oar , the defined corotation 
radius Rcr and its error estimate ARcRioi each galaxy. Also given are parameter 1Z, error 
estimate ATZ and Hubble stage T (from HyperLeda). 
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is not well-suited for modelling the dust lanes, but indeed 
we got a better fit for these features with fast-rotating bar. 
However, the spiral response with such pattern speed choice 
is quite bad. In the case of NGC 4995 the fit is not very good 
and a value giving a fast rotating bar cannot be excluded. 
NGC 6384 has very long spiral arms, which are completely 
reproduced with the listed value of 1Z. Only the innermost 
part can be reproduced with a fast rotating bar. 

An alternative explanation to slow bars, that the spiral 
structure in late type galaxies typically rotates with a lower 
pattern speed than the bar, and that our method gives too 
much emphasis to the spiral, is certainly possible for some 
galaxies of the sample (NGC 3726, NGC 4303, NGC 4995, 
NGC 6384). However, there are several galaxies (NGC 289, 
NGC 578, NGC 1832), where also the fit to the innermost 
features favours low pattern speed. If the bar in these galax- 
ies is fast rotating, then even the innermost part of the spiral 



rotates with a lower pattern speed, and the bar may be rel- 
atively insignificant in their dynamics. 



4.3 Size and strength of the bar 

The comparison of the determined corotation radius to 
the galaxy size, the latter estimated to be proportional to 
isophotal radius R25, does not show any clear trend with 
galaxy morphology. Regardless of the Hubble type, Rcr can 
be anything between 0.3 and 1.0 times 7?25- The situation is 
different if we compare the dimensionless pattern speed TZ 
to relative bar size, measured by the estimated semi major 
axes Rbar, and scaled to the isophotal radius R25 (Fig. [S] 
top. Now there is a clear correlation: the slow bars tend to 
be also the smallest ones: with one exception, the galaxies 
with Rbar/ R25 ^ 0.25 have 1Z ^ 1.4. When all the galaxies 
are plotted in Rcr/R25 vs. Rb ar /R25 one can see for the 
smallest bars (Rbar/ R25 « 0.15 — 0.35), the corotation reso- 
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Table 2. Comparison sample. 
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Figure 7. Top: a comparison sample of previous determinations 
of bar pattern speeds as a function of 1Z. Filled circles show pat- 
tern speed estimates based on Tremaine- Weinberg method, tri- 
angles estimates based on modelling, squares estimates based on 
Puerari-Dottori method and diamonds estimates based on mor- 
phological or kinematical arguments. Bottom: The results of our 
sample compared with results obtained with Tremaine- Weinberg 
method. 



Figure 8. Top: best fitting 1Z vs. relative size of the 
bar Rbar/^25, bottom: relative corotation radius RcbJ R25 
VS.R bar /R25- 

bars. There is no clear correlation between pattern speed 
and spiral strength. 



nance is close to 0.4 it^s- This corresponds to t ypica l values 
for the radii, where lElmegreen fe Elmegreenl (Il995l ) found 
the end of the symmetric two-armed spiral to be in non- 
barred galaxies and galaxies with small bars. For the larger 
bars, the bar size essentially correlates with the corotation 
radius (Fig. [8] botto m), but the scatter is large. 

In some studies l|Patsis et al.lll99ll ; iKranz et all 120031 ) 
it has been suggested that in non-barred galaxies the end 
of the spiral or at least the end of its symmetric part lies 
near corotation or inner 4/1-resonance. The possible inde- 
pendence of corotation resonance radius from the bar size in 
the galaxies with smallest bars could suggest that the spiral 
in these galaxies indeed has a different pattern speed, whose 
corotation radius is regulated by the galaxy size. 

I n Fig. [9] w e comp are 1Z to parameters Qb and Qs 
from iButa et al.l (|2005l ). which characterize the strengths 
of bar and spiral components (maximum tangential force 
divided by the azimuthally averaged radial force in each ra- 
dius). With one exception, galaxies with slow bars (1Z > 
1.4) have Qb ^ 0.3. However, there are several galaxies 
with similar bar strength but fast bar. The combined re- 
sults with bar size and stren gth are in accordance with 
lElmegreen fc Elmegreenl (|l995l ) who found by morphologi- 
cal arguments that the flat (and typically long) bars tend to 
be faster rotators than the exponential (and typically short) 



4.4 Bulge-to-total flux ratio, galaxy luminosity 
and color 

The top frame of Fig. \W\ shows R as a function of 
bulge-to-total fl ux fraction (B/T-ratio) in H-band (see also 
ISalo et alj 2007). The apparent bulge contribution is one of 
the criteria in Hubble classification, so it is not surprising 
that we get almost a mirror image of Fig. [6] 

We studied the dependency between TZ and absolute B- 
luminosity (c alcula t ed fro m the B-magnitudes in RC3) and 
distances by iTullvl (|l98ct ). but found no clear correlation. 
The same was the sit uation with absolu te if-magnitudes 
based on 2MASS data l|Jarrett et al.|[2003l ) and neither there 
was correlation with colour. However, if we plot the galaxies 
in (B/T, A')-plane f Fig. 1101 bottom), then the fast and slow 
galaxies roughly occupy different areas: slower bar favouring 
lower bulge fraction and/or higher galaxy luminosity. 

4.5 Pattern speed vs. dust lane morphology 

Accor ding to gas dynamical simulations by lAthanassoulal 
(1992), the existence of leading offset dust lanes in the 
bars tells us two things: that there is an ILR and that 
1Z — 1.2 ± 0.2. The first condition comes from the existence 
of so-called X2-orbits, which require the presence of an ILR: 
the dust lanes form as the orientation of the orbits of gas 
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Figure 9. Top: best fitting 1Z vs. bar strength Qb, bottom: 1Z 
vs. spiral strength Qs- 



Figure 10. Correlation with bulge flux fraction. Top: 7Z vs. B/T, 
bottom: fast (black circles), slow (open circles) and intermediate 
(grey circles) bars in (B/T, ii")-plane. 



clouds change from parallel orientation of zi-orbits (with 
respect to the bar) in the outer parts of the bar to the per- 
pendicular orientation of a;2-orbits in the inner parts. The 
second condition comes from the gas dynamical simulation, 
only in this parameter region the shape of the shock corre- 
sponds to the observed dust morphology. 

We have checked the dust lane morphology of the sam- 
ple galaxies from OSUBSGS images and images from various 
other sources, e.g. NED and the web page of de Vaucouleurs 
Atlas of Galaxies (jButa et al.ll2007h . Omitting cases where 
the inclination makes recognition and classification of dust 
lanes ambiguous, we found that more than third of the galax- 
ies in the sample does not show classic dust lane morphology. 
In some of the cases the dust lanes are missing altogether, 
perhaps due to lack of interstellar matter in the bar region, in 
others the morphology is more complicated or even chaotic. 
Of the seven galaxies with lZ m in > 1.4, the dust lane mor- 
phology is uncertain due to inclination for NGC 3726 and 
NGC 4995. In NGC 289 the dust lane morphology is not of 
the "classic" type, in NGC 578 it is more or less chaotic. 
NGC 1832 and NGC 6384 do not show dust lanes in the bar 
region. In NGC 4303 the dust lane morphology is of classic 
type suggesting that the bar in this galaxies may be fast 
rotating. Indeed, if we compare a model for NGC 4303 with 
1Z ~ 1 (Fig. [3]), then the gas morphology in the bar region 
shows features that look a bit like the observed dust lanes, 
but the spiral structure is not well reproduced. 



5 CONCLUSIONS 

We have modelled 38 barred galaxies with simulations us- 
ing collisionless test particles and inelastatically colliding 
test particles. The gravitational potentials used were derived 
from //-band images from Ohio State University Bright Spi- 
ral Galaxy Survey. 

Our findings are as follows: 

1) the average value of the dimensionless pattern speed, 
1Z = Rcn/Rbar depends on the morphological type, being 
about 1.2 for the early type barred galaxies (SBO/a-SBab), 
1.4 to intermediate type (SBb) and 1.8 to late type (SBbc- 
SBc). 

2) When error estimates are considered, all early type 
spiral galaxies in the sample (T ^ 2) are in accordance with 
having fast bars (1Z < 1.4). 

3) If the derived pattern speed corresponds to that of 
the bar, then intermediate and late type galaxies include 
both fast and slow bars. In five cases the bar remains in 
slow bar domain even if a "worst case scenario" allowed by 
the error estimates in pattern speed and bar size determina- 
tions is considered. However, the existence of multiple pat- 
tern speeds, i.e. spiral rotating more slowly than the bar, is 
a viable alternative in some cases, especially for NGC 4303. 

4) Slow bars are short when compared to isophotal ra- 
dius 7?25- Especially, with only one exception, galaxies with 
Rbar/R25 ^ 0.25 have a slow bar. 

5) With one exception, the slow bars have bar strength 
Qb ^ 0.3, although several bars with similar strength are 
also found to be fast. 
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6) There is no clear correlation with the absolute B- or 
if-magnitudes of the galaxies or the B — K color index. 

7) Galaxies with fast and slow bars seem to occupy dif- 
ferent regions when plotted in a two-dimensional coordinate 
system defined by bulge-to-total flux ratio and absolute K- 
magnitude. 

8) Omitting possible systematic errors, the error es- 
timates of the model-based pattern speeds are typically 
smaller than with Tremaine- Weinberg method. If the results 
by these two methods are found to be consistent, then mod- 
elling would enable pattern speed estimation for large galaxy 
samples. 
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